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Abstract 

The physics of interacting nuclear spins arranged in a crystalline lattice is typically 
described using a thermodynamic framework [1 J : a variety of experimental studies 
in bulk solid-state systems have proven the concept of a spin temperature to be not 
only correct [2, 3j but also vital for the understanding of experimental observations |4j. 
Using demagnetization experiments we demonstrate that the mesoscopic nuclear spin 
ensemble of a quantum dot (QD) can in general not be described by a spin tempera- 
ture. We associate the observed deviations from a thermal spin state with the presence 
of strong quadrupolar interactions within the QD that cause significant anharmonic- 
ity in the spectrum of the nuclear spins. Strain-induced, inhomogeneous quadrupolar 
shifts also lead to a complete suppression of angular momentum exchange between 
the nuclear spin ensemble and its environment, resulting in nuclear spin relaxation 
times exceeding an hour. Remarkably, the position dependent axes of quadrupolar 
interactions render magnetic field sweeps inherently non-adiabatic, thereby causing 
an irreversible loss of nuclear spin polarization. 
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The study of nuclear spin physics by optical orientation experiments in bulk semi- 
conductor materials has been an active field of research over the last decades [3J El E] ■ These 
research efforts have shown that using the electron as a mediator, it is possible to transfer an- 
gular momentum from light onto nuclei, thereby establishing a nuclear spin polarization that 
is orders of magnitude higher than the equilibrium nuclear polarization at cryogenic tem- 
peratures. As a result, the effective nuclear spin temperature in such an optically pumped 
system can be pushed to the low mK regime. Combining these optical pumping schemes 
with nuclear adiabatic demagnetization techniques borrowed from bulk NMR experiments [3] 
would be a natural extension to these experiments that could lead to a significant further 
reduction of the nuclear spin temperature. This approach, previously demonstrated in bulk 
semiconductors [51 E], suffers from the fact that in most systems where optical orientation 
of nuclear spins is possible, nuclear spin relaxation is too fast to allow for a significant re- 
duction of magnetic fields in an adiabatic way. Here, we use the exceedingly long nuclear 
spin relaxation time in self-assembled QDs [D] to implement an "adiabatic" demagnetization 
experiment on the system of ~ 10 5 nuclear spins. 

The mesoscopic ensemble of nuclear spins in a QD can be conveniently polarized and 
measured by optical means [5J El [101 El H3] • To this end, we use the photoluminescence (PL) 
of the negatively charged exciton (X -1 ) under resonant excitation of an excited QD state. 
It has been shown previously [H] that under appropriate excitation conditions, 20 — 50% of 
the QD nuclear spins can be efficiently polarized in a timescale of a few ms. The resulting 
dynamical nuclear spin polarization (DNSP) can then be measured through a change in the 
Zeeman splitting, AEos, of the X -1 recombination line [H]; this energy shift due to the spin 
polarized nuclei is commonly referred to as the Overhauser shift (OS) [13] . 

A remarkable feature of the QD nuclear spin system is the excellent isolation from its 
environment if the QD is uncharged. Fig.[1Ji shows the corresponding free evolution of 
Pnnc (proportional to Ai^os) i n a QD subject to an external magnetic field B ext = 2 T. 
The nuclear spin relaxation time clearly exceeds one hour and does not vary appreciably 
over the magnetic field range relevant to this work [9] . Since the bulk material surrounding 
the QD remains unpolarized during the experiment (see Methods), the long nuclear spin 
lifetime indicates that nuclear spin diffusion between the QD and its environment is strongly 
suppressed. We attribute this quenching of spin diffusion to the structural and chemical 
mismatch between the InGaAs QD and its GaAs sourrounding [T31 US] • The very slow 
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nuclear spin relaxation leaves room for further manipulation of the QD nuclear spin system 
after optical pumping. In particular, we can study how P nuc behaves under slow variations 
of external parameters and thereby study the validity of spin thermodynamics for the QD 
nuclear spin system. 

If the QD nuclei were describable using a thermodynamic approach, P nuc would be 
aligned with the external magnetic field _B ex t and would be described by Curie's law 
7-Pnuc = B cxt C/T spin [3] (here, 7 is the nuclear gyromagnetic ratio, C the Curie constant 
and T spin the nuclear spin temperature). An adiabatic lowering of B ext from an initial value 
Bi to a final value B{, would conserve P nuc and lead to a reduction of T spin by a factor 
Bf/B v In general, cooling by adiabatic demagnetization is possible for any system where 
the spin entropy S is conserved and a function of -B ex t/T spm only. The ultimate limit to 
the achievable cooling is determined by nuclear spin interactions which give the dominant 
contribution to S at low magnetic fields. The strength and nature of these interactions can 
be phenomenologically described by a random local magnetic field B\ oc . In most cases, B\ oc 
is given by the nuclear dipolar couplings (~ 0.1 mT). As soon as B ext « B\ OC) the local fields 
randomize an established nuclear spin polarization and thereby limit the efficiency of the 
adiabatic spin cooling to B\ oc /B{. The resulting behavior of nuclear spin temperature and 
polarization as a function of Bf is sketched in Fig.[T]D: for B cxt = 0, the spin temperature 
remains finite and the nuclear spins are completely depolarized. Amazingly, this depolar- 
ization is a reversible process, provided that S is a conserved quantity at all fields. When 
the spins are re-magnetized to a magnetic field exceeding B\ oc , their polarization recovers 
along the direction of the magnetic field and in particular conserves the sign of its initial 
spin temperature. 

To test the validity of spin thermodynamics for the QD nuclear spins and to study the 
possibility of adiabatic cooling in this system, we performed demagnetization experiments 
on a QD, as illustrated in Fig.[]J. A circularly polarized "pump" pulse of length r pump is 
used to polarize the nuclear spins. We then linearly ramp B cxt from B\ to B{ with a rate 
7b = 10 niT/s. At the final field Bf, the remaining degree of nuclear spin polarization is 
measured using a linearly polarized "probe" pulse of length r pro b e [9] . This experiment is 
repeated at various values of Bf to record the process of "adiabatic" (de)magnetization. 

Figure[TJi shows the result of a demagnetization experiment performed on the nuclear spin 
system of an individual QD. The nuclei are polarized with a pump pulse r pump = 300 ms 
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at Bi = 1 T and measured at B{ with a probe pulse r pro b e = 5 ms. At a rough glance, 
this measurement qualitatively follows the behavior depicted in Fig.[T]D. A closer inspection 
however, reveals significant deviations: upon ramping the external field to Bf = — 1 T we 
only recover 63% of the initial P nU c- In addition, by measuring P nuc (Bf) we determined the 
value of the local field to be B\ oc = 290 mT: this value is about three orders of magnitude 
larger than typical nuclear dipolar fields. Finally, we observe that even for Bf = 0, the QD 
has a remnant nuclear spin polarization P^™. To verify that we do not induce an unwanted 
increase of spin entropy by sweeping B ext too fast, we repeated our experiment for values of 
7b of 5,2.5 and 0.8 mT/s (crosses in Fig.JTji). Within the experimentally accessible range, 
7b has no influence on our observations. 

The discrepancy between our experimental findings and the predictions from a thermo- 
dynamical treatment of nuclear spins becomes even more pronounced if we increase P nU c(-Bi) 
(which can be achieved by increasing B\ [T2l [T3]). Fig. [2^, shows an experiment where we 
demagnetize the polarized nuclear spins starting from B x = 2.2 T to a final field B' f = — 1 T 
(black data points). We then reverse the sweep direction of the magnetic field and ramp 
_B ext back to Bi (gray data points). This experiment shows a considerable hysteresis of 
the nuclear spin polarization as a function of Bf. In particular, -P^u™ changes sign for the 
two sweep directions of B ext . Furthermore, the magnitude of P^ 1 , resp. the width of the 
observed hysteresis curve depends linearly on the initial degree of nuclear spin polarization 
and on B x (Fig.ga). 

To obtain more information about the source of irreversibility of P nuc during magnetic 
field sweeps, we performed a further experiment, where we optically orient the nuclear spins 
at B\ = 1 T, ramp the field to a value B\ < Bi and then back to B\ = Bf where we measure 
the remaining degree of nuclear spin polarization. The result of this experiment (Fig.^) 
indicates that the magnetic field sweeps start to induce irreversibilities in P nuc as soon as 
|Pcxt| £ Pioc « 300 mT. 

Finally we note that the experimental observations described here do not depend on the 
sign of the initial nuclear spin temperature (T sp i nj i). We have repeated the demagnetization 
experiments for T sp i n; i < (i.e. cr~ laser excitation at B\ > 0, not shown here) and observed 
values of -Bi oc and P^ consistent with the measurements presented in Fig.[T] and Fig.[2j 
These measurements are complicated by the fact that for T spin < 0, nuclear spin pumping 
is rather inefficient [13] , leading to a low degree of DNSP and therefore a smaller signal to 
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noise ratio than for T spin > 0. 

The three principal features of our experiments, the existence of -P^™, the hysteretic 
behavior of P nuc and the partial irreversibility of our demagnetization experiment, result 
from a violation of the nuclear (Zeeman) spin temperature approximation [HE]. We explain 
these features by taking into account the strong inhomogeneous quadrupolar interactions 
(QI) of the nuclear spins in a QD [T7J [HI [19]. The self-assembled growth of InGaAs QDs 
is driven by a strong lattice-mismatch between InGaAs and its surrounding GaAs matrix, 
which results in a heavily strained QD lattice. As a consequence, QD nuclei experience 
large electric field gradients which couple to the nuclear quadrupolar moment. The resulting 
quadrupolar Hamiltonian [20] , 

#Q = ^(^-3 /(/ + 1)) ' (1) 

is characterized by a nuclear quadrupolar frequency vq (proportional to the local strain at 
the nuclear site) and a quadrupolar axis z' (with corresponding unit vector e z > along the 
main axis of the local electric field gradient tensor), i is the nuclear spin angular momentum 
operator with quantum number I and I z i = I • e z i . For typical strain values of 2 % |22j , we 
find Uq ~ 2.8 MHz for As and 1.2 MHz for In[23J. For comparison of the interaction 
strength of Hq with a pure nuclear Zeeman Hamiltonian H% = 7I • B ext , it is convenient 
to express the QI strength by an equivalent magnetic field Bq = huq/^. For As and In, 
we find Bq = 388 mT and 125 mT, respectively; the corresponding mean value agrees well 
with our experimental estimate for B\ oc . 

The spectrum of a nuclear spin with quadrupolar frequency vq depends strongly on 
the angle 9 between e z > and the external magnetic field (directed along e 2 ). Figure^ 
shows the Eigenenergies of a nuclear spin with I = 3/2, as a function of B ext /BQ. At 
-Bext = 0, the spectrum is governed by Hq, which pairs the nuclear spin states into doublets 
with angular momentum projections ±m z > on e 2 /. The doublets are split by an energy 
\huJm, ,,m t+i\ — ( m z' + 1/2)/^q, respectively. Conversely, in a high magnetic field, the 
spectrum is determined by H% with nuclear angular momentum being quantized along the 
axis e z . Even at arbitrarily high fields however, the spectrum is significantly perturbed by 
.f/Q and never becomes perfectly harmonic. 

We modeled our demagnetization experiment using the steady state solution of a rate 
equation for the populations p\ m \ of spin states \m), which are mutually coupled through 
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dipolar interactions (Fig.[3p and c). The nuclear spins are initialized with a Boltzmann 
distribution at B ext = (see Methods) and the evolution of the p\ m \s is calculated as a 
function of B ext . Due to the unequal nuclear spin level spacings, only nuclear spin flip- 
flops that preserve p\ m \ (V|m)) are energetically allowed in general and therefore the spin 
populations remain invariant as a function of B ext . Varying B ext will change the relative 
nuclear spin level spacings in the nonlinear way depicted in Fig.[3j3. Since the P| m )'s remain 
invariant as B ext is reduced, the nuclear spins are driven into a state which is out of thermal 
equilibrium (i.e. not Boltzmann-distributed) . At specific values of B ext (red markers in 
Fig.^), transition energies between distinct pairs of nuclear spin states can coincide — a 
situation denoted as a "cross-over" of nuclear spin transitions [lj. At those fields, the p| m )'s 
are no longer constant and the nuclear spin levels involved in the cross-over can relax to 
a Boltzmann distribution. The irreversibility observed in our magnetic field sweeps is a 
consequence of this partial relaxation of nuclear spins to thermal equilibrium. We speculate 
that the resulting increase of the nuclear spin entropy is induced by an energy-conserving 
coupling to the environment of the nuclear spins. If the minimal energy gap of the anti- 
crossing induced by the dipolar coupling between two interacting nuclear spins at their 
cross-over is smaller than the coupling to the environment, pure dephasing of the nuclear 
spin transitions will induce irreversible cross-over transitions and S will increase. 

Upon sweeping B ext through zero (red box in Fig.^), dipolar interactions will couple 
the states m z i = ±1/2. The associated passage through the avoided crossing between these 
single-spin states is adiabatic and preserves the respective populations in the two lowest- 
lying spin states. In contrast, nuclear dipolar interaction can not couple any of the states 
with \m z i\ > 1/2 due to conservation of energy and angular momentum. The spin states 
in the \m z /\ = 3/2-manifold will therefore cross and in particular preserve their populations 
P3/2 and P-3/2- The imbalance between these populations (p 3 / 2 < P-3/2 in Fig.[3| will result 
in a remnant polarization -P^™, even if B ext is strictly zero. 

We averaged our model over a set of parameters 6 and uq to account for the strong 
inhomogeneity of QI over the QD (see Methods). The result of this full simulation, is 
shown in Fig.[3]i. We highlight that the good qualitative agreement with our experimental 
results (Fig. [2^l) is rather insensitive to the set of parameters used in our simulation. In 
particular, the choice of the distribution for the parameters 9 and uq did not affect our 
results significantly. Furthermore, our simulation treats the QD spin system as a pure spin- 
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3/2 system, while for In, / = 9/2. A numerical treatment of the full InGaAs nuclear spin 
system is beyond the scope of this paper and would most likely not alter the qualitative 
behavior of our simulations (see Methods). 

Our results show that the nuclear spin system of a self-assembled QD provides a rare 
example for a solid-state nuclear spin ensemble that can not be described by a nuclear spin 
temperature [21] . We note that if one could assign a spin temperature to the QD nuclear 
spin system, optical pumping combined with adiabatic demagnetization of the nuclear spins 
would be a novel and efficient way of nuclear spin cooling in QDs without QI: possible systems 
include nuclear spin- 1/2 systems, such as 13 C-nanotube QDs [22], where QI is inherently 
absent, or strain-free semiconductor nanostructures [26], such as epitaxially grown droplet- 
QDs [27J . There, adiabatic nuclear spin cooling would only be limited by nuclear dipolar 
interactions resulting in I?i oc « 0.1 mT. Achieving nuclear spin cooling to temperatures 
« 100 nK should be feasible in these systems, opening ways to studying the remnants of 
nuclear magnetic phase transitions in the mesoscopic system of QD nuclear spins |28j . 
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Individual QDs were studied using the photoluminescence (PL) of X~ x under resonant 
excitation of an excited QD state. The QD sample was grown by molecular beam epitaxy 
on a (100) semi-insulating GaAs substrate. The approximate composition of the QDs after 
self-assembled growth and post-growth annealing was Ino.5Gao.5As. For individual optical 
addressing, the QDs were grown at a low density of < 0.1 fim~ 2 . The QDs were spaced by 
25 nm of GaAs from a doped n ++ -GaAs layer, followed by 30 nm of GaAs and 29 periods of 
an AlAs/GaAs (2/2 nm) barrier which was capped by 4 nm of GaAs. A bias voltage applied 
between the top Schottky and back Ohmic contacts controls the charging state of the QD. 
Optical pumping of QD nuclear spins was was performed at the center of the A -1 stability 
plateau in gate voltage, where PL counts as well as the resulting OS were maximized [14"]. 

The QD sample was immersed in a liquid Helium bath cryostat equipped with a super- 
conducting magnet and was held at the cryostat base temperature of 1.7 K. The PL emitted 
by the QD was analyzed in a 750 mm monochromator allowing for the determination of 
spectral shifts of the QD emission lines with a precision of ~ 1 /xeV [12] • A combination 
of an optical "pump-probe" technique, together with linear ramps of the applied magnetic 
field were used to adiabatically demagnetize the QD nuclear spins (see Fig.JT^); technical 
details of the pump-probe setup are given elsewhere [9] . The "pump" pulse consists of a 
circularly polarized laser pulse of duration r pump which is used to optically orient the QD 
nuclear spins [13]. We typically achieve an OS of AEqs = 60 /xeV at B t — 1 T, corresponding 
to nuclear spin polarization P nuc « 35% or 7] ~ 1.5 mK (for £>; = 2.2 T, AE os = 89.5 fieV 
and P nuc ~ 50%). In the range of B ext relevant to our experiment, P nuc oc B x [T3] such that 
the initial nuclear spin temperature 7] is roughly constant and on the order of few mK [6] 
for all values of B x . 

Directly after applying the pump pulse to the QD, the gate voltage is switched to a value 
where the QD is charge-neutral. In this regime, nuclear spin polarization has an exceedingly 
long relaxation time on the order of hours [9] (see Fig.[T^). We note that we can exclude any 
significant nuclear polarization of the bulk material surrounding the QD. The observation 
of DNSP in our experiment depends sensitively on the excitation laser energy, which we 
tune to an intra-dot (p-shell) excitation resonance with a width of ~ 300 fieV and located 
« 36 meV above the PL emission energy. The sharpness and energy of this excitation 
resonance makes any excitation processes which involve the creation of free electrons in 
the bulk very unlikely [2T] . Furthermore, the pumping time T pump = 600 ms used in our 
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experiment is much too short to lead to a significant bulk nuclear spin polarization, even if 
some free electrons were created during laser illumination. 

Details of the model 

The model we developed to explain our experimental findings is based on the steady state 
solution of a rate equation for the populations of a nuclear spin 1 = 3/2 system. The nuclear 
spins are initialized with a Boltzmann distribution over the spin states at B ext = B\. The 
assumption of a thermal distribution of nuclear spin levels at B ext = B x is justified by the 
fact that nuclear spins are polarized by hyperfine interaction with the QD electron: optical 
pumping of the electron leads to a broadening of its spin states by several fieV [32], allowing 
for electron-nuclear flip-flops between the electron and any two given nuclear spin states 
which are coupled by the hyperfine interaction. It is therefore reasonable to assume that 
the occupations of nuclear spin levels at B\ follow a Boltzmann distribution. 

We then change the magnetic field by keeping the populations of spin levels fixed. Only 
at the specific fields where cross-relaxation is permitted (Fig.J^), we allow for a local ther- 
mal equilibrium to be established between the spin levels involved in the cross-relaxation 
transitions. All other populations and the total energy of the nuclear spin system remain 
constant. Upon sweeping through B ext = 0, we assume that the levels m z i = ±1/2 undergo 
an adiabatic passage through an anticrossing induced by the coupling of these two states by 
dipolar interactions. Spin states with m z > = ±3/2 however remain uncoupled and undergo 
an adiabatic level crossing which preserves their populations. 

The result of our simulations is shown in Fig.|3]3 and d of the main text. We illustrate 
the evolution of the occupations of the individual nuclear spin states in Fig. [3^, where we 
show the spectrum of a nuclear spin for the parameters z/q = 3 MHz, 9 = 0.37r/2 and 
7 = 10 MHz/T. The occupations of the individual levels is encoded by the thickness and 
gray shade of the corresponding lines. Magnetic fields where cross-relaxation processes 
take place are indicated by red lines. We repeated this calculation for a set of angles 6 £ 
{0.1, 0.2, 0.6}| and quadrupolar frequencies z/q G {—4,-3,-2,2,3,4} MHz over which 
we average our results. Since the local strain in our QDs can be both tensile and compressive, 
positive and negative values for z/q are possible. By solving the complete Hamiltonian 
Hq + Hi, we can relate the occupancies of the spin levels to our experimentally observed 
nuclear spin polarization — the expectation value (I z ) of the nuclear spin polarization along 
the direction of -B e xt- Fig.[3jl of the main paper shows the result of our simulation in form 
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of the calculated evolution of P mic function of B{. 

We note that our model is a great simplification of the actual experimental situa- 
tion. First, we completely neglect cross-relaxation events between nuclei of different 
(6, vq) — values. Second, our calculation was performed for a spin-3/2 system for simplicity, 
while the actual QD nuclear spin system consists of a mixture of spin-3/2 (Ga, As) and spin- 
9/2 (In), which further complicates the situation. While a numerical treatment of the full 
InGaAs nuclear spin system is beyond the scope of this paper, we argue that such a treat- 
ment would not alter the physical picture conveyed by our simulation. Including 1=9/2 spins 
would lead to a nuclear spin spectrum similar to the one illustrated in Fig.|3]3. The number 
of magnetic field values where cross-relaxation events would be energetically allowed would 
increase compared to the case of 1=3/2, but these events would still be singular in the sense 
that for most values of B ext , the nuclear spins could not thermalize. The system would thus 
still be driven out of thermal equilibrium and the relaxation events during cross-relaxation 
would lead to an increase of nuclear spin entropy. Including flip-flop events between In and 
As nuclear spins would have a similar effect: these transitions would be allowed for a subset 
of close nuclei and would allow for partial thermalization only at specific values of B ext . 
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FIG. 1: Demagnetization of QD nuclear spins, (a) Free decay of P n uc at -B ex t = 2 T for 
an uncharged QD after optical pumping of the nuclear spins for r pump = 600 ms. The gray curve 
shows an exponential decay with a time constant of 1 h for comparison, (b) Theoretical prediction 
of nuclear spin temperature and polarization during adiabatic demagnetization from a field (red 
arrow) to Bf. (c) Schematic of the experimental procedure for adiabatic demagnetization of QD 
nuclear spins. The nuclei are optically pumped at B ext = B{ (T sp i n .i ~ mK). Directly after the 
pumping pulse, the electron is ejected from the QD. B eyLt is then linearly ramped at a rate 7b to 
a value Bf at which we measure P n uc- (d) Experimental (de)magnetization of QD nuclear spins. 
B[ = 1 T as indicated by the red arrow, 7# = 10 mT/s and AEos(B{) = 57 fieV. The gray curve 
is a fit according to the theoretic predictions shown in (b); we find B\ oc = 290 mT. Blue, green and 
red crosses show a similar experiment, with 75 = 5, 2.5 and 0.8 mT/s, respectively {B^ = 0.5 T 
for these data-points). 
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FIG. 2: Irreversibility and hysteresis in the demagnetization experiment. Black circles: 
Same experiment as in Fig.[T]i, with = 2.2 T as indicated by the red arrow (AEos(B{) = 
89.5 fieV). After reaching B' { = —1 T, we reverse the magnetic field sweep direction and bring 
the nuclei back to the initial field (gray circles), (b) The remnant nuclear spin polarization P^utf 
(normalized to the value P nuc (B[) found in (a)) as a function of B\. As P nuc (Bi) oc B[, the nuclear 
spin temperature after optical pumping is roughly constant for all values of Bi in this measurement, 
(c) After polarization of nuclear spins at B^ = 1 T (red arrow), we sweep B cxt to £?[ and then back 
to B[, where P n uc is measured. The magnetic field sweeps become partly irreversible as soon as 
|-Bf| < 0.3 T« Bq. The lines in the figures are guides to the eye. 
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FIG. 3: Modeling of the demagnetization experiment. Local electric field gradients induced 
by strain in self-assembled QDs result in strong QI for the nuclear spins, (a) Model of local strain 
axis-distribution e z > within a QD. (b) Spectrum of nuclear spins (1=3/2) under the influence of 
both Hq and Hz for a variety of angles between e z / and e z . (c) Simulation of QD nuclear spin 
demagnetization for a particular setting v q = 3 MHz and = 0.157T. Nuclear spin populations p 
are represented both by line-thickness and grayscale of the lines that indicate the energy of the 
nuclear spin states. At Bi = 1 T, the nuclei are initialized with a Boltzmann distribution over 
their spin states. The populations remain constant for most values of Bf. Only if a cross-over 
of nuclear spin transitions occurs (red markers for B{ > 0), the occupations of the involved spin 
states evolve to a (local) thermal equilibrium distribution (see text). We simulate this process for 
a set of configurations {9, vq} and calculate the corresponding magnetization P nuc oc (I z ). (d) 
shows the resulting nuclear spin polarization as a function of B{ starting at B{ (red arrow), which 
qualitatively reproduces the experimental findings shown in Fig. [2] 
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